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Abstract 

We show that for certain configurations of two chemostats fed in parallel, the presence of two different 
species in each tank can improve the yield of the whole process, compared to the same configuration 
having the same species in each volume. This leads to a (so-called) "transgressive over-yielding" due to 
spatialization. 
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1 Introduction 

The study of biodiversity and environment on the functioning of ecosystems is one of the most important 
issues in ecology nowadays (see for instance [H]). 

Recent investigations on the functioning of microbial ecosystems has revealed that combinations of certain 
species and resources enhance the functioning of the overall ecosystem, leading to a "transgressive over- 
yielding" ( |13| [8]). In those experiments, the environment is modified by altering the composition of the 
resources and its richness ([16). In the literature in ecology, the consideration of spatial heterogeneity in 
the environment (niches, islands...) is often promoted as a leading factor for the explanation and prediction 
of abundance and distribution of species and resources in ecosystems. Comparatively, the influence of the 
spatialization on the performances of the functioning of ecosystems has been much more rarely studied ([I])- 
From the view point of the industry of biotechnology and waste water treatment, the spatial decoupling 
of bio-transformation is known to impact significantly the yield conversion, and engineers are looking for 
design of bio-processes that offer the best performances (see [HI [TJ ID ) ■ For instance, a series of continuous 
stirred bioreactors is known to be potentially more efficient than a single one (for precise distributions of 
the volumes) in [6J [5j [10] . Recent investigations on the chemostat model have revealed that other simple 
spatial motifs, such as parallel interconnections of tanks, can also present interesting benefits ([3]). In this 
framework, most of the studies are conducted for a single species. Only a few of them investigate the role of 
biodiversity and spatial heterogeneity on the performances of the ecosystem (|15[ ITTj). 

In the present work, we show that spatialization can be another factor responsible for a transgressive 
over-yielding. For this purpose, we consider the usual framework of the chemostat and show that splitting 
one volume into two tanks in parallel with different dilution rates can enhance the performances at steady 
state when two different species are present in each vessel, compared to the situations of having the same 
species present in both volumes. Here, optimal performance is achieved when the value of the nutrient at 
steady state in the outflow of the system is minimized. This notion is inspired by the scenario of waste 
water treatment processes where the input nutrient to the reactor is the pollutant, and the purpose of the 
treatment plant is to minimize the pollutant concentration. 
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We first study the case of a single species and show that two reactors in parallel instead of one does 
not bring any benefit, although bypassing one tank may improve the performance in some cases. Then, 
we give sufficient conditions for the growth curves of two species that guarantee the existence of two tanks 
configurations that offer a transgressive over-yielding. 

2 One species 
2.1 One reactor 

Let a chemostat be operated with constant volumetric flow rate Q, constant volume V, and constant input 
nutrient concentration S m . We call D = Q/V the dilution rate. In case only 1 species is present, the 
equations for the nutrient concentration S(t) and the species concentration x(t), are: 

S = D{S in -S)-^x (1) 
7 

x = x(f(S)-D) (2) 

where 7 G (0, 1] is the yield of the conversion of nutrient into species biomass, and f(S) is the per capita 
growth rate, a smooth, increasing function (with f'(S) > for all S > 0) with /(0) = 0. We also assume 
that / is concave, i.e. f"(S) < for all S > 0. Finally, we let D* = f(S in ). 

A typical example satisfying all these conditions would be the Monod function f(S) — mS/(a + S) for 
positive a and m. 

The first statement of the following result is well-known (see for instance the textbook [33]). 
Proposition 1 Let x(0) > 0. Every solution (S(t),x(t)) 0/ (jT|) — (|2j) satisfies: 

lim S(t) = X(D), 

t— >+oo 

where 

X(D) = f-\D), if D < D*, and =S m , if D > D* (3) 

The (continuous) function X(D) is increasing on [0,D*) (since A' > there) and nondecreasing on [0, +00). 
X(D) is convex if D 6 [0,D*) (since X" > there), but not if D £ [D*,oo). 

Proof . That X(D) is continuous, is obvious from the definition, and that it is increasing on [0, D*) follows 
from the fact that it is the inverse function of the increasing function f(S), S £ [0, S m ). To show the 
convexity of X(D) for D e [0,D*) we implicitly differentiate = S, S E [0,5 m ], twice: 

(r l )V(s))f(s) = 1 

(/- 1 )"(/(5))(/'(5)) 2 + (/- 1 )'(/(5))/"(5) = 0, 

and thus that 

ff—iyi _ (/ ) / 
U ) - (/ , )2 

Since /' > and hence also {f^ 1 )' > 0, and since /" < 0, it follows that (J -1 )" > 0, hence / _1 is convex, 
and therefore so is X(D) as long as D e [0, D*). Nonconvexity of X(D) when D G [0, +00) follows by picking 
an arbitrary D > D* , and letting D* = pO + (1 — p)D for an appropriate p S (0, 1). It follows that: 

X{ P + (1 - p)D) = X(D*) = S m 

> (1 - p)S m = (1 - p)D = P A(0) + (1 - p)X(D), 
and thus X(D) is not convex if D G [0, +00). | 
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2.2 Two reactors 



Assume that we split the inflow channel to the chemostat in two channels, and the reactor volume in two 
parts. The first reactor is fed by the first channel with volumetric flow rate aQ for some a e [0, 1] and it 
has volume rV for some re [0,1]. The second reactor is fed by the second channel with volumetric flow rate 
(1 — a)Q, and it has volume (1 —r)V. Note that the total reactor volume of both reactors equals the initial 
single reactor volume. Note also that the initial inflow channel has simply been split into two channels. Of 
course, we assume that for both species into bioreactors optimal environmental conditions (temperature, pH 
etc.) arc applied. 

To determine the asymptotic value of the nutrient concentration in each of the two reactors, we use the 
Proposition from the previous section. This value is completely determined by the dilution rate that each 
reactor experiences. For the first reactor, this dilution rate is 

aQ _ a D 
rV r 

whereas it is 

(l-q)Q = 1-a 
(l-r)V l-r 
for the second reactor. Therefore, the asymptotic nutrient values are 




in the first and second reactor respectively. Assume that we mix the outflows of both reactors (in the same 
proportions a and 1 — a as the volumetric flow rates at their respective input channels). The asymptotic 
value of the nutrient in the mixture is then given by the function 

F(a,r):=aA(^D)+(l-a)A(^£>), (4) 

where (a, r) e (0, l) 2 . Clearly, this function is smooth in (0, l) 2 . We would like to extend it to [0, l] 2 . First 
we extend it continuously to (a, r) ^ (0, 0) and ^ (1, 1) by defining 

F(a,0) := aS m + (1 - a)A ((1 - a)D) for a > 0, 
F(a,l) := a\{aD) + (l-a)S m for a < 1 

Finally, we continuously extend it to [0, l] 2 by defining: 

F(0,0)=F(1,1) ~X(D) 

An important property of F is that: 

F(a 7 r) = X(D), whenever a = r. 

This is not surprising since all configurations with a — r result in the same dilution rate D for each of the 
two reactors, which is the same as the dilution rate of the initial single reactor. The asymptotic nutrient 
concentrations in each reactor will be the same, and equal to X(D), which equals the asymptotic nutrient 
concentration of the single reactor. The mixing of these concentrations with proportions a and 1 — a, yields 
X(D). In other words, for these configurations, the splitting of the reactor is merely artificial, and the 
asymptotic outcome is the same as if the reactor would not have been split at all. 

An obvious question is to determine those configurations (a, r) where F(a, r) is minimal, and also to see 
if these minima are lower than the corresponding value X(D) which is obtained in case of a single reactor 
(or in case of any configuration satisfying a — r in view of the remark above) 
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We first define a convex polygonal region in the unit square [0, l] 2 : 

C = {(a, r) e [0, l] 2 | (a/r)D < D* 
and ((l-a)/(l-r)) £><£>*} 

and we note that C is a nontrivial set (i.e. has nonzero Lebesgue measure) if and only if 

D < D*, 

an assumption we make in the remainder of this paper. This condition on D simply expresses that in case of 
a single reactor, the asymptotic nutrient concentration X(D) will be strictly less than the maximal possible 
value S m . Next, we also define two triangular regions in [0, l] 2 : 

Ti = {(a, r) G [0, l] 2 | (a/r)D > D*} 

T 2 = {(a,r)6[0,l] 2 | ((1 - a)/(l - r)) D > D*)} 

and we note that C, T\ and T 2 are pairwise disjoint, and that they cover [0, l] 2 . 
The following Lemma is usefull for the following. 



6 : (a,r) e [0,1] x (0,1] ^ — (5) 

r 



Lemma 1 The function 

is convex on its domain. 
Proof . One has 

and furthermore 
It follows that 

Saa > , S rr > , 5 aa S rr - 8 2 ar = . 
Thus the Hessian of 8 is positive semi-definite. | 
We introduce the following function 

T in (D) = A(.D) + D\'(D) (6) 

that plays an important role in the following. 

The following Proposition shows that the performance of a single species chemostat can only be improved 
if S m is small enough, and this is achieved by bypassing the second reactor with a unique, specific flow rate. 

Proposition 2 The restriction of F to the convex set C, F\c, is convex. Moreover there holds that 



min F = 

[0,1] = 



X(D) - F(a,a), Va E [0,1], for S m > T in (D) 
min Q F(a,0) < X(D), for S in < T in (D) 



Proof . Let us first show that F\c is convex. We consider the functions L and R defined on C: 



L(a, r) = aX (?-D) , R(a, r) = (1 - a)X {^-^ D 



(7) 



and show that L and R are convex on C. 
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Let (ai,ri), (a 2 ,r 2 ) belong to C and /j be a number in [0, 1]. One has 
T = fiL(ai,ri) + (1 - (i)L(a 2 ,r 2 ) 

= [i ai \ (j^d) + (1 - M)"2A (j^d) 

(JLOL\ 



= (fiai + (1 - (i)a 2 ) 

By convexity of the function A, it follows 

T > + (1 - fi)oe 2 )\ 

= {mi + (! - m)«2)a 



/xai + (1 — ^t)a2 \ ri 



D 



{l-(l)a 2 W^?_£> 



/Ltai + (1 - jU)d!2 \ r 2 

/xaf/ri + (1 - ^)allr 2 D 

\1<X\ + (1 - (l)OL 2 

fj,S(ai,ri) + (1 - fj)8(a 2 , r 2 ) 



D 



(ia\ + (1 - (i)a 2 

The function A being increasing and 6 being convex (see the former Lemma), one obtains 



r > ((mxi + (1 - (i)a 2 ). 

( 8(nai + (1 - (i)a 2l [iri + (1 - (x)r 2 ) 



A 



D 



[lax + (1 - (i)a 2 

= [pai + (l-(t)a 2 )X( -7- r — TJ> 

V m + (1 - M)^2 J 

= L(/xai + (1 - fi)a 2 , (in + (1 - (i)r 2 ) 

Similarily, it can be shown that R is convex on C. The function F being the sum of the two functions L and 
R, one concludes that F is also convex on C. 



For (a, r) G int(C), we have that 



aD 



A' ( —D 



1-a 
1-r 



-D 



(l-a)D (1-a 



1-r 



A 



1 



F r = - 



a 2 D x , ( ol 
— r — A 



-D 



(1 - a) 2 T> (1-a 
A 



1-r 



(8) 
(9) 
(10) 



(1-r) 2 

We note from © and (fit))) that every point on the diagonal where a = r in int(C) is a critical point, i.e.: 

^a\a—r — -^r\a—r 0; 

and thus, since F is convex in C, it follows that 

mmF = A(-D) = F(a,a), Va e [0, 1] (11) 

Since [0,1] 2 is the disjoint union of C, T\ and T 2 , the proof of the Proposition will be complete if we 
study the minimisation of F on Ti and T 2 . Notice that T 2 = {(a,r) G [0, l] 2 | (1 - a, 1 - r) G Ti} and 
F(a, r) = F(l — a, 1 — r). Consequently, one has min^ F = minT 2 F, and we can study the restriction of F 
to Ti only. When (a, r) G Ti, then 



F(a,r) = aS m + (1 - a) A 



1-a 

1-r 



-D 
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Note that F has no critical points in int(Ti) because 



(1 - r) 2 \l-r 

and thus the extrema of are necessarily located on the boundary of T\. We calculate the values of F 
on the various parts of the boundary of T± , and bound them below: 

F(a,0) = aS in + (1 - a)A((l - a)D) 
F(l,r) = S m >X{D) 

F(a,a^) = aS m + (1 - a)A ( - \~ n\ 



D* / 

> aS in + (1 - a)A((l - a)D) = F(a, 0) 

where we used that A(.) is non-decreasing. One has 

F a (a,0) = S in - A((l - -DA'((1 -a)D)(l -a) 

^aa(a,0) = 2DA'((l-a)L>)+D 2 A"((l-a)D)(l-a) 

and one concludes that a i— > F(a, 0) is strictly convex, A being a convex increasing function on the domain 
[0, D]. Furhermorc, one has 

F a (0,Q) = S in -T m (D) 
F a (l,0) = S m 

Consequently the unique minimum of F(a, 0) over a is reached at a = exactly when F a (0, 0) > 0, that is 
to write S m > T m (D). Finally, notice that one has F(Q, 0) = X(D) and the proof of the Proposition now 
follows from (fTT|) . | 



Remark 1 When Si n < T m (D) the best configuration is obtained for the by-pass of one tank with a different 
than or 1 . For the Monod function, the threshold on S m has the following expression 

„s F)a ( m \ 

T m (D) = 1+ . 

m — D \ m — D J 



3 Two species 

We consider two species with respective per capita growth rate functions f\ and fi, satisfying the conditions 
of a single species we imposed in the previous section. We assume that there is a unique S € (0, S m ) such 
that 



fi(S) > f 2 (S)i£SG(0,S) 
f 2 (S) > f^S) if 5 6 (5, +oo) 

and we let 

D = MS) = f 2 (S), D\ = A(S in ), D* 2 =f 2 (S m ) 

We assume as before that 

D < mm(Dl,D* 2 ) 



G 



Definition 1 We say that a configuration (a,r) £ (0, l) 2 corresponds to transgressive overyielding if 

(12) 



G(a,r) := aAi (^Dj + (1 - a)A 2 (jz% D> \ 
< mm(F 1 (a,r),F 2 (a,r)), 



where Xi is defined as in ([3]) but with f = fi and D* = D* , and Fi is defined as in dU) but with A = Xi . 

Let us first consider the case where a = r. Then 

G(a,a) = aXi{D) + (1 - a)X 2 (D) 
Fi(a,a) = aXi(D) + (1 - a)Xi(D) = Ai(D) 
F a (a, a) = aX 2 (D) + (1 - a)X 2 (D) = X 2 {D) 

Consequently, since any convex combination of two non-negative numbers cannot be smaller than the 
smallest of these two numbers, it follows that 

G(a,a) = aXi(D) + (1 - a)X 2 (D) 

> mm(Xi(D),X 2 (D)) =mm(Fi(a,a),F 2 (a,a)), [ ' 

which implies that no configuration with a = r can correspond to transgressive overyielding. Moreover, if 
Ai(D) X 2 (D), then the inequality in (|13|) is strict for all a £ (0,1). By continuity of G,Fi and F 2 this 
inequality remains strict for all configurations (a,r) near such configurations where a — r and a s (0, 1), 
implying that we will not find configurations corresponding to transgressive overyielding nearby. On the 
other hand, if Xi(D) = X 2 (D) = S, or equivalently if D = D, then the inequality in (|T51) reduces to an 
equality for all a £ (0, 1). In this case we show that there always exist configurations nearby corresponding 
to transgressive overyielding. 

The two next Propositions show that there always exist configurations of two parallel chemostats with 
only one of the species in each reactor volume that perform better than in the case in which for the same 
configuration, only one of the species is present in both reactors, regardless of the selected species. 

Proposition 3 Let D = D and fix a £ (0, 1). Then there always exist configurations (a,r) near (a, a) that 
correspond to transgressive overyielding. 

Proof . We start with the simple observation that for any (a, f) £ (0, l) 2 : 

— — 1 ;. 1 ~ ^ > i 
f > \ _ f < 

Therefore, for fixed a £ (0, 1), we can always find (a,f) £ (0, l) 2 near (a, a) such that 

<3 ^ ^ 1 — ot - 
-D < D < -D 

r 1 — f 

Moreover, since D = D < min(_D^, D 2 ), we may assume w.l.o.g. that 

1 - a - 

D < min(D*DZ) 

I _ f 

The properties of f\ and f 2 imply that: 

< Ai (jD^j < X 2 (jDj < S 
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and thus in particular that: 

< Ai (~D\ < A 2 (jD 

. ( 1 — a -\ , / 1 — a - 

< A2 {— D ) < Al ( — D 

There holds that 

G(a, f) = SAi (jDj + (1 - a)X 2 

Fi(5,f) = aAi (jD\ + (1 - 5)Ai 

F 2 (a,f) = &\ 2 (jD^j + (1 - a)X 2 (j^D 

Noting that each of G(a, f), F\ (a, f) and F 2 (a, f) is a convex combination with the same weights a and 1 — a 
of a pair of distinct positive numbers from a collection of four positive numbers that are ordered as indicated 
by (33]), it follows that: 

G{a,r) < mm(F 1 (a,r),F 2 (a,r)) , 
and thus the configuration (a,r) corresponds to transgressive overyielding. | 

The case where D ^ D, can be handled similarly: 
Proposition 4 Let D ^ D. Then there exist (a,r) € (0, l) 2 such that 

G. - 1 — OL 

-D < D < D < mm(D*D*) (14) 

r 1 — r 

and this configuration (a, r) corresponds to transgressive overyielding. 

Proof . Using the fact that D,D < min(-D|, -D|), it is not hard to show that there always exist (a,r) € 
(0, l) 2 such that (14J) holds. Indeed, the half spaces {(a,r)\a/r < D/D} and {(a,r)|(l — a )/(l — r) < 
min(£>*, D^/D} are easily seen to intersect in the open unit square (0, l) 2 . From this follows that (TT4|) 
holds with (a,f) replaced by (a, r), and then the rest of the proof of the previous Proposition can be 
carried out with (a, f) replaced by (a, r). In conclusion, the configuration (a, r) corresponds to transgressive 
overyielding. | 

4 Numerical examples 

We consider two species whose graphs of growth functions /j (i = 1, 2) intersect away from zero : 

/lO) = T— ' ^ 2 ( s ) = n q' S , 

3 + s 0.3 + s 

One can easily check that fi(S*) — f 2 {S*) for 5* = 2.4. For the simulations, we have chosen D = 0.9 that 
is slightly above fi(S*) = f 2 (S*) — 0.88 (see Figured]). For this value of the dilution rate D, the break-even 
concentrations of the two species are 

A X (D) = 2.499, A 2 (D) = 2.993 

The thresholds given by formula ([6]) are 

T{ n (D) = 7.08 , T 2 in {D) = 35.84 
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Figure 1: Growth curves and dilution rate D 
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40 


10 


5 


mm a Fx (at, 0) 


2.50 (= 


= Xx(D)) 


2.50 (= Xx(D)) 


2.35 


min Q F 2 (a, 0) 


2.99 (= 


= X 2 (D)) 


2.21 


1.61 



Table 1: Best values of by- pass configurations. 



We have studied numerically three values of S m , larger than both thresholds, in between or lower. The best 
values for the bypass configurations are reported on the following table For S m = 40, the best configuration 
is obtained for a = 0.41 and r = 0.36 that achieves the minimum of G equal to 0.234, which represents a 
gain of 12% compared to having the best species (the first one) in a single tank. Accordingly to the Table 
[TJ the best configurations for S m = 10 and S m — 5 are obtained for the bypass of the first tank, using 
the second species only in the second tank (although this species has the worst break-even concentration). 
Nevertheless, we show on Figures [2] and [3] the benefit of the parallel configuration with non-null volumes for 
different distributions of the volumes given by the parameter r. 

One can check that in any case the minimum of G(-,r) over a is lower than min(Ai(Z?), \2(D)) and is 
reached for a value of a close to r. 



5 Conclusion 

Our main contribution is to show that whatever are monotonic growh rate functions of two species, whose 
graphs have an intersection away from the origin, there always exist configurations of two parallel chemostats 
with only one of the species in each reactor volume that perform better than in the case in which for the 
same configuration, only one of the species is present in both reactors, regardless of the selected species. We 
beleive that this result is of interest for potential applications in bio-industry. Further investigations could 
be conducted to determine more precisely what gain could be expected. Our simulations show that this gain 
could be significant. 
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Figure 2: Graphs of Fi(-,r) and G(-,r) for S in = 10 




Figure 3: Graphs of Fi(; r) and G(; r) for S ln = 5 
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